Fringe spacing and phase of interfering matter waves 
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We experimentally investigate the outcoupling of atoms from Bose-Einstein condensates using 
two radio- frequency (rf) fields in the presence of gravity. We show that the fringe separation in 
the resulting interference pattern derives entirely from the energy difference between the two rf 
fields and not the gravitational potential difference. We subsequently demonstrate how the phase 
and polarisation of the rf radiation directly control the phase of the matter wave interference and 
provide a semi-classical interpretation of the results. 
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Ever since the first realisations of Bose-Einstein con- 
densates (BECs) in dilute atomic gases, their coherence 
properties have been the subject of much investigation. 
The first clear demonstration that BECs possess long 
range phase coherence was through the interference of 
two spatially separated condensates [l| . Since then, other 
experiments have studied the coherence properties us- 
ing atom laser output from an array of tunnel coupled 
condensates in an optical standing wave 2] , Bragg spec- 
troscopy P(, density fluctuations 3 and intereferometry 
An elegant scheme to probe condensate coherence, 
based on interfering atom laser beams, was reported by 
Bloch et al. (3,0- This used two radio frequency (rf) 
fields to outcouple atoms from different locations within 
a condensate. A high contrast matter wave interference 
pattern was observed at temperatures well below the 
BEC transition temperature, confirming the phase co- 
herence of the condensate. A numerical model of two 
outcoupled modes agreed with the experimental obser- 
vations 0. More recently, the atom- by-atom build up 
of a matter wave interference pattern has been observed 
using single atom detection [T(|. To date, experiments 
have primarily focussed on the visibility of the interfer- 
ence patterns. In this paper, we describe experiments 
which address the fringe spacing, phase and nature of 
the interference. 

Outcoupling atoms from Bose-Einstein condensates 
with rf fields has been used extensively to produce 
beams of atoms, generally referred to as "atom lasers" 
[Til l~L2L H^ . The rf radiation drives resonant (stimu- 
lated) transitions from a trapped Zeeman sublevel to an 
untrapped state in which the atom falls under gravity. 
Outcoupling occurs at locations where the total energy 
difference between the trapped and untrapped states is 
equal to Huj r f. This is usually determined by the Zeeman 
potential so that the resonant condition may be written 

hu; r f = iib9f\B(t)\ (1) 

where \ib is the Bohr magneton, gp is the Lande ^-factor 
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and B(y) is the magnetic field. 

Consider a condensate trapped in a cigar-shaped mag- 
netic potential of the form U (r) = muo 2 (n 2 x 2 + y 2 + z 2 )/2 
where m is the mass of the atom, uo z = uo y is the trap- 
ping frequency in the tight directions of the trap and 
k = uo x /uo z . Atoms can be outcoupled from the surface 
of an ellipsoid of the magnetic equipotential which sat- 
isfies the resonance condition QJ. However, gravity will 
cause a displacement of the minimum of the total po- 
tential from the magnetic field minimum. This gravita- 
tional sag means a harmonically trapped condensate will 
be displaced from the magnetic field minimum by a dis- 
tance, zo = —g/uo 2 , where oj z is the trapping frequency 
in the direction of gravity. This displacement is typically 
greater than the size of the condensate, so that the el- 
lipsoidal equipotential surfaces can be approximated by 
planes which intersect the condensate at different heights, 
z. In this situation, the dependence on the x and y coor- 
dinates can be neglected for many quantitative purposes 
and only the z dimension need be considered. 

In previous work 0,0? two rf fields of frequencies, uj\ 
and 6^2, were used to outcouple atoms from a condensate. 
The two spatially separated resonances were interpreted 
as creating two slits from which atoms were extracted 
from the condensate. The outcoupled atoms formed two 
matter waves which interfered, in close analogy with a 
Young's double slit experiment. The visibility of the in- 
terference pattern provided a measure of the first order 
phase coherence of the condensate. 

The outcoupling points, z\ and z^ used in 0, 0, 
were chosen to be centred around the middle of the con- 
densate, located at Zq, the minimum of the combined 
magnetic (harmonic) and gravitational (linear) potential. 
Under this condition, the gravitational energy difference 
between the two outcoupling points, determined by the 
slit separation, Az = z\ — 22, is exactly equal to the dif- 
ference in energy between the two applied rf fields. This 
can easily be seen from the derivative of the magnetic po- 
tential, where AE « muo 2 zAz. At the central position, 
zo, we find AE = H(ui — UJ2) = mgAz. 

However, this result is only true when z = zo, (where 
z = (z\ + Z2) /2, is the distance from the magnetic 
field minimum to the centre of the two slits). If the 
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FIG. 1: (a) Total energy (solid lines) of atoms in the trapped, 
ttif — —1 and untrapped ttif — states, for the parame- 
ters used in our experiments. The dashed and dash-dotted 
lines represent two pairs of rf fields, with equal Acu, used to 
outcouple atoms from the BEC at different locations. The 
gravitational energy difference mgAz may vary for a fixed 
Acj, however, when the interaction energy (shaded) of atoms 
in the untrapped state is included, the total energy difference 
between the two outcoupled matter waves is always equal to 
HAlj. (b) The measured outcoupled matter wave beams dis- 
play an interference pattern with a constant fringe spacing 
X(z) given by equation 0] for a fixed Au (in this case equal to 
2tt x 1000 s _1 ), independent of Az. 

two resonant points are not symmetrically located about 
the centre of the trap, the z 2 dependence of the mag- 
netic potential means that the slit separation for a fixed 
Auo = uj\ — U02 varies inversely with z. Thus the gravita- 
tional energy difference between the two resonance posi- 
tions, is not necessarily equal to HAuj and Az may change 
significantly across the width of a condensate. For a har- 
monic potential, the slit separation is approximately 



In a Young's double slit experiment, the fringe spacing, 
A, of the interference pattern is proportional to Az -1 . 
However, in dual rf out coupling experiments with a fixed 
Acj, A is independent of Az. While the gravitational 
energy difference, mgAz, between the two resonant loca- 
tions can change, this is not the only energy to consider. 
The fringe spacing of the interference pattern depends 
on the total energy difference between the two indistin- 
guishable outcoupling paths and must always equal HAuj 
to satisfy energy conservation. This is a general result 
which we discuss below, for a BEC in the Thomas-Fermi 
(TF) regime. 

Consider the specific case of an F = 1 8T Rb TF con- 
densate. The total energy of the trapped state \F = 
1, = —1) consists of the sum of its magnetic, gravita- 
tional and mean field energies. The untrapped (outcou- 



pled) state 1 1,0) experiences negligible magnetic poten- 
tial, but, while still within the condensate, experiences 
both the mean field and gravitational potentials. In ID, 
the total energy of a particular substate can be written 

as 

E mF (z) = -m F {]^mu 2 z z 2 + ii B 9fBq) - mgz + #i d |^0)| 2 

(3) 

where Bq is the magnetic field at the minimum of the 
trap, gid is the ID effective interaction strength (as- 
sumed to be the same for all m^, which is approximately 
true but not an essential point in this discussion) and 
1^(^)1 2 = ^2m F \ ri l ) m F { z )\ 2 is the- total atomic density. 
The energies of the trapped |1,— 1) state and the un- 
trapped 1 1,0) state are plotted (solid lines) in Fig. da) 
for the parameters used in our experiments. The shaded 
regions indicate the mean field contribution to the total 
energy. Also shown are two pairs of rf fields, dashed and 
dash-dotted lines, with the same Acj, chosen to lie within 
the width of the condensate. As the two pairs are cen- 
tred around different z, the resulting Az for each pair is 
different, but the total energy difference is HAuj. 

In the TF limit the interaction energy between the con- 
densate and the outcoupled state exactly compensates 
for the difference in gravitational potential at different 
slit locations. The density profile, \i(j(z)\ 2 , mirrors the 
shape of the magnetic trapping potential so that the en- 
ergy splitting between the two states is always given by 
the difference in their magnetic potentials. Additionally, 
the energy of trapped atoms within the condensate is 
independent of z so that only the final energies on the 
ttif = curve determine the energy difference between 
the two outcoupled beams. 

Having established that the fringe spacing, A, depends 
only on Acj, it can easily be shown that 

m = (4) 

where Af = Aco/2tt. A varies with z because the outcou- 
pled atoms accelerate in the z-direction under gravity, as 
can be seen in Fig. IHb). 

We have performed a range of experiments to verify 
this for several values of Aco with the resonant points 
centred around various z positions. Our experimental 
procedure for producing condensates has been described 
elsewhere 14] and was used here with only slight mod- 
ifications. An atom chip is used to produce near pure 
condensates containing 2 x 10 5 87 Rb atoms in the |1, — 1) 
ground state. Our chip design facilitates the production 
of relatively large condensates in highly stable trapping 
fields. The final trapping frequencies are 160 Hz in the 
tight direction and 6.7 Hz in the weak direction. The 
elongated geometry of the trap means we must cool well 
below the 3D critical temperature T c to produce fully 
phase coherent condensates (typically for our param- 
eters is less than T c /2 |15j). Outcoupling is induced by 
turning on two rf fields of the same amplitude, with fre- 
quencies lji and UO2 tuned to be resonant with atoms in 
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the condensate, and Rabi frequencies, Q = HB9F\B\/2h, 
of 50 Hz for each rf source. After outcoupling for 10 ms, 
the trap is left on for a further 3 ms before being turned 
off abruptly. An absorption image is taken after 5.3 ms 
of free expansion. 

We first checked the reproducibility of the fringe spac- 
ing for a fixed Aco of 2ir x 1000 s -1 as Az was varied 
over the range 390 nm to 560 nm by varying z. While the 
visibility of the interference pattern decreased near the 
edges of the condensate, the wavelength was consistent 
with the value predicted by equation (|3J to within 1%. 
Experimental values of A were determined by taking an 
image of the out coupled atoms, similar to that shown in 
Fig. IHb), converting the z spatial axis into a time axis 
through the relation t = \Jl(z — zo)/g, integrating the 
output over x and fitting a cos 2 function to the data. 
The uncertainty in A is determined by the uncertainty in 
the fitted frequency of the cos 2 function. Variations at 
the level of 1% are within our experimental uncertainties 
and not significant when compared to what would be ex- 
pected if A was determined by Az -1 . This would lead to 
a variation of more than 30% across the range of values 
we measured. 

A semi-classical interpretation of these dual rf out- 
coupling experiments, based on the interference of the 
applied rf fields, may also be used to understand these 
experiments. Defining the mean frequency, Co = (uo\ + 
6^2 )/2, and the beat frequency, 5 = (uj\ — c^2)/2, and 
recalling the standard trigonometric identity 

smujit + sina^ = 2smOt cos St (5) 

we see that the sum of two oscillating fields is equivalent 
to an amplitude modulated carrier wave. In our case, the 
carrier frequency, u), is typically three orders of magni- 
tude higher than the beat frequency, S. Adding a phase, 
</>, to one of the rf fields shifts the phase of both the 
carrier and beating terms by half this amount. 

We may now consider the condensate interacting with 
a single rf field at the carrier frequency, which is ampli- 
tude modulated in time. The number of atoms outcou- 
pled is proportional to the Rabi frequency squared (ie. 
proportional to the amplitude squared of the rf field at 
time t) and is modulated at 25 = Auo. Once outcoupled, 
the atoms fall under gravity and, provided they have a 
low spread of initial momenta, the outcoupled density 
will be modulated in time. 

To demonstrate this, a sequence of dual rf outcoupling 
experiments was performed using fixed values uo\ and 002 
(Auo = 2tt x 500 s -1 ) but varying the relative phase of 
the two rf fields. This has the effect of shifting the phase 
of the rf beat note. All other experimental parameters 
were kept fixed. Two examples of the data obtained are 
shown in Fig. [21 Absorption images of the outcoupled 
atoms appear on the left and the beat note of the corre- 
sponding rf fields used for outcoupling (measured on an 
oscilloscope) are shown on the right. The z axis of the 
absorption images has been rescaled by \[%z — zo)/g to 
linearise the time axis for ease of comparison with the rf. 




FIG. 2: The phase of the interfering matter wave beams is 
determined by the phase of the beating rf fields used in to 
drive the outcoupling. (a) and (b) represent different runs of 
the experiment under identical conditions apart from differ- 
ent phases of the applied rf field. On the left are absorption 
images of the condensate (top) and outcoupled atoms and 
on the right is the beat note of the corresponding rf used to 
drive the outcoupling, measured on an oscilloscope. The ver- 
tical axis indicates the time before the image was taken and 
was obtained for the atom images through the relationship, 
t = ^2{z - z )/g 



It is clear from these images that the outcoupled atoms 
correspond to the largest amplitude of the rf beat note. 

In order to quantify this, we have analysed a series of 
similar data in which the relative phase of the rf fields 
was allowed to vary randomly over the range to 2tt. 
The fitting procedure described earlier was applied to 
all of the absorption images to determine the phase of 
the outcoupled beam. A similar fit was applied to the 
square of the measured rf beat note and the two phases 
are plotted against each other (filled sircles) in Fig. 03 
The dashed line through this data is a plot of y — x. The 
phase of the modulated atom beam matches very well the 
phase of the applied rf field. 

In these experiments both rf fields were provided by 
passing the two rf currents through the same coil. This 
means the rf field in the vicinity of the BEC was lin- 
early polarised, perpendicular to the quantisation axis x. 
We have also performed experiments using separate, or- 
thogonally mounted coils where each rf current was sent 
through a different coil. The rf interference is no longer 
linearly polarised but rather a field whose polarisation 
varies from vertical linear, to left hand circular, to hori- 
zontal linear, to right hand circular and back to vertical 
linear in a single beat period (T = l/Aw). This is anal- 
ogous to the optical field used in lin _L lin sub-Doppler 
(Sisyphus) laser cooling but the field is periodic in 
time rather than space. 

The outcoupling transition from |1,— 1) to |1,0) re- 
quires H-fi of angular momentum, and can only be driven 
by <j + radiation. When the rf is derived from a single 
coil, it is easy to see that the maximum amplitude of 
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FIG. 3: Plot of the fitted phase of the modulated atom beam 
against the phase of the beating rf field. Filled circles rep- 
resent experimentally measured data using a single coil and 
the dashed line is a plot of y = x. Open squares are exper- 
imentally measured phases when the rf is produced by two 
separate (near) orthogonally mounted coils and the dotted 
line is a plot of y = x + 0.557T. 

cr + radiation occurs when the two rf fields are in phase. 
The linearly polarised field may be decomposed into two 
counter rotating circular fields and it is the <r + compo- 
nent of this which couples to the atom. With perpendic- 
ularly oriented coils however, the maximum amplitude of 
the cr + field occurs during the circular polarised phase of 
the beat note. This happens when the rf fields are 7r/2 
out of phase. 

Also shown in Fig. [3| (open squares) is a plot of the 
phase of the atom laser versus the phase of the beating 
rf fields performed with (near) perpendicularly oriented 



coils. The atomic output is phase shifted by approx- 
imately 7r/2 from the rf as expected. The slight mis- 
match between the measured shift of 0.557T and the ex- 
pected shift of 0.57T was due to imperfect alignment of 
the rf coils (precise perpendicular alignment would have 
impeded optical access in our setup). For coils mounted 
antiparallel maximum outcoupling would occur when the 
two rf fields are tt out of phase. 

In conclusion, we have studied the origins of the fringe 
spacing and phase of matter wave interference patterns, 
produced by outcoupling atoms from a BEC with two 
rf fields. We have shown that the energy difference be- 
tween the two rf fields determines the spacing of the in- 
terference pattern, not the gravitational potential differ- 
ence determined from the classical slit separation. Semi- 
classical arguments based on interfering rf fields correctly 
predict the experimental observations. These also show 
how the phase and polarisation of the rf field determine 
the phase of the observed matter wave interference pat- 
tern. This extends previous work which looked at the 
fringe visibility for the specific case where HAcu = mgAz 
0. Our findings do not contradict the phase coherence 
studies reported in ref. [jj. Indeed, any random phase 
gradients within the BEC would lead to random initial 
velocities that would degrade the observed interference 
patterns. Finally, we note that we have also performed 
experiments with cold thermal atoms and see (as in Q) 
that the visibility of interference pattern diminishes, due 
to the thermal spread of velocities in the trapped gas. 
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